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An unpolarized charge current passing through a chaotic quantum dot with spin-orbit couphng 
can produce a spin polarized exit current without magnetic fields or ferromagnets. We use random 
matrix theory to estimate the typical spin polarization as a function of the number of channels in 
each lead in the limit of large spin-orbit coupling. We find rms spin polarizations up to 45% with 
one input channel and two output channels. Finite temperature and dephasing both suppress the 
effect, and we include dephasing effects using a new variation of the third lead model. If there is only 
one channel in the output lead, no spin polarization can be produced, but we show that dephasing 
lifts this restriction. 
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I. INTRODUCTION 

The generation and control of spin polarized currents, 

\ in particular without magnetic fields or ferromagnets, 
is a major focus of recent experimental and theoreti- 
cal work. This includes the spin Hall effect, which pro- 
duces spin currents transverse to an electric field in a 
two-dimensional electron system (2DES) with spin-orbit 
coupling, with spin accumulation at the edges. ^ Simi- 
larly, the magnetoelectric eflFect^'^'^ produces a steady 

\ state spin accumulation when an electric field is applied 
to a 2DES with spin-orbit coupling. The accumulation 

' can be uniform^i^ in the case of uniform Rashba spin- 
orbit coupling^ or at the edges of a channel in either the 
Rashba model^i^ii^ or with spin-orbit coupling induced 
by lateral confinement.— Experiments have observed 
current induced spin polarization in n-type 3D samplesi^ 

\ and in 2D hole systemsiiii^ with spin polarization esti- 
mated to be up to 10%."'^^ Further work suggests a spin 
polarized current can be produced by a quantum point 
contact (QPC) with spin-orbit coupling fi^i- in a car- 
bon nanotubefi^ in a ballistic ratchet, in a torsional 
oscillator^ or in vertical transport through a quantum 

■ well,^! 

Here we show that generating a polarized current from 
an unpolarized current is a generic property of scattering 
through a mesoscopic system with spin-orbit coupling. 
We propose using many-electron quantum dots (outside 
the Coulomb blockade regime) with spin-orbit coupling 
to produce partially spin polarized currents without mag- 
' netic fields or ferromagnets. Due to the complicated 
boundary conditions of the quantum dot, we do not solve 
for the spin polarization in terms of any particular spin- 
orbit coupling model, geometry, and contact configura- 
tion. We estimate the effect for a ballistic system in the 
limit of strong spin-orbit coupling by performing a ran- 
dom matrix theory (RMT) calculation for the spin polar- 
ization, allowing consideration of realistic quantum dot 
devices robust to details of shape and contact placement. 
Finely tuned systems should be able to exceed these po- 
larizations, but these results provide a useful benchmark 
for whether a particular tuned system is better than a 




FIG. 1: A tuned model quantum dot with A'^ = 1 channels 
in the left lead and M = 2 channels in the right lead. A 
skew scatterer sends all z spins to the top channel and all 
— 2 spins to the bottom channel. The shaded area in the 
bottom channel has Rashba spin-orbit coupling of precisely 
the strength to rotate a down spin at the Fermi energy to 
an up spin, thus producing a spin polarized exit current from 
any input current, while respecting time reversal symmetry. 

generic chaotic one. We use a density matrix formalism 
throughout, which allows us to develop straightforwardly 
a spin-conserving dephasing probe, using a new variant 
of the third-lead technique for accounting for dephasing. 
Dephasing and finite temperature both reduce the ex- 
pected polarization. Without dephasing, we find that if 
there is only one outgoing channel then no spin polariza- 
tion is possible, which was first shown by Zhai and Xui^ 
Interestingly, with dephasing, spin polarization can be 
produced with only one outgoing channel. The case of 
polarized input currents will be discussed elsewhere^SS 

Analogous calculations have been performed by Bar- 
darson, Adagideli, and Jacquod in a four-terminal geom- 
etry, to study the transverse spin current produced by an 
applied charge current.— 



II. SETUP AND SYMMETRY RESTRICTIONS 

We consider non-interacting electrons in a quantum 
dot with two attached leads connected to large reservoirs. 
For any electron current entering from the leads, we can 
describe the output state in the leads in terms of the S- 
matrix of the dot, including any tunnel barriers between 
the leads and the dot. We assume negligible spin-orbit 
coupling in the leads and consider the lead on the left 
(right) to have N (M) spin-degenerate channels at least 
partially open at the Fermi energy, and let K = N + 
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M. As usual, the channel wavefunctions are normalized 
so all channels have the same flux. The S-matrix S is 
a 2K X 2K unitary matrix of complex numbers. For 
spin 1/2 particles with spin-orbit coupling, however, it 
is convenient to consider S to he a K x K matrix of 
quaternions. We give a brief introduction to quaternions 
and explain why they are easier to work with. 

We choose a representation of the quaternions such 
that a quaternion g is a 2 x 2 matrix of complex num- 
bers q = g*^°H2 + « ^''^■'o'a" where are the Pauh 
matrices and g^''-' € C. We define three conjugates of q: 
complex conjugate q* ~ g("^*l2+i J2 quaternion 
dual q^ = q^^^t2 — i^l'^^'^'^tn and Hermitian conjugate 
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For a. K X K matrix of quaternions, Q, we de- 
fine {Q% = Q,;, (Q«)^^. = Q/, and (Qt) = {Q* f. 
By convention, the trace of the quaternion matrix Q is 

Given a K x K quaternion matrix Q, we can associate 
it with a 2K x 2K complex matrix A in the obvious 
way. Note then that is equivalent to A^, the usual 
Hermitian conjugate of a complex matrix, but Q* is not 
equivalent to A*. Note further that the trace convention 
implies that tiA = 2trQ. 

The quaternion representation is convenient, as the 
time reversal operation for a scattering matrix can sim- 
ply be written as S" ^ S^.^^ The S-matrix of a system 
with time reversal symmetry (TRS) is self-dual. 



If w" 



[w 



is the K X K quaternion density ma- 



trix of the incoming (outgoing) current. 



The density matrix describing the unpolarized incoherent 
combination of all N incoming channels is 



2N I Om 



(1) 



That is, w'" = Pl/2N where Pl is the projection onto 
the channels of the left lead. We choose trw"^ = 1/2, due 
to the trace convention. 

The Landauer-Biittiker formula gives the conductance 
in terms of the S-matrix. We write the K x K quater- 
nion S-matrix as (j */) with r (r') being the N x N 
(M X M) reflection matrix and t {t') the M x N {N x M) 
transmission matrix. Then we write the Landauer- 
Biittiker formula in units of 2e^ jh as 



G = tr(tt^), 

= 2NIy{PrSw'''S'^), 
= 2iVtr(Pfli(;™*), 



where Pr is the projection onto the channels of the right 
lead. Since is normalized to represent one input par- 
ticle entering the system, g = 2tr(PflZ(;°"*) is the proba- 
bility for that particle to exit through the right lead. The 
conductance is TV times this probability, so we call g the 
conductance per channel in the left lead. 



Similarly, we define a vector spin conductance22i (i.e., 
exit spin current divided by voltage) in units of e/27r 
as 



= 2Nix{aPRW°''^). 



Then 



g'' = 2tr(CTPKW°"' 



(3) 



(4) 



is the spin conductance per channel in the left lead. 
Hence, 5* is the /z-component of the spin polarization 
of the exit current times the probability of exiting into 
the right lead. Thus, the spin polarization of the current 
in the right lead is p = <7^/<7j with \p\ < 1. 

We can, of course, construct g, g^, andp using only the 
S-matrix and not the density matrices w"^ and it;°"' . The 
density matrix approach, however, gives the flexibility to 
consider arbitrarily correlated states of incoming current 
and also to look for arbitrary correlations in the out- 
going current. ^■^ We will also use it to straightforwardly 
derive a method of accounting for non-magnetic dephas- 
ing in a device with spin-orbit coupling. To complete the 
translation to the standard notation of conductances, we 
consider sending up- or down-polarized electrons into a 
sample and collecting either up- or down-polarized elec- 
trons, giving a conductance matrix^^ 



Q ^ / G|| G|x 
\Gn Gil 



(5) 



with the total charge conductance being G = G^ | +G|| + 
G|| 4- G||. Gct.o-' is the conductance for an input current 
of spin a' and an exit current of spin cr, for cr, a' =t, i- We 
translate the quaternion representation into the standard 
notation by noting that the up-polarized incoherent input 
current has input density matrix = ^^^Pl- The 
output density matrix is 711°"' — Sw^-^S'^ and the portion 
representing the output in the right lead is t ^^^^ . The 
Landauer-Biittiker formula gives, in units of 2e^ /h, 



Gtt =iVtr(Pfli^w^^* 



=tr( 



l + g.-! j- 1 + 0-3 



tt). 



(2) Similarly, 



Gil 
G 
G 



1-0-3 f l-t-oa f\\ 
2 ^ 2 h 

2 



tr(-^t 

Ti =tr( 2 . 2 



Xi =tr(i^il^tt) 



(6) 



(7) 
(8) 
(9) 



from which we see that G = tr(tt^^), which is the usual 
Landauer-Biittiker formulai^ 

Though there are several spin-orbit coupled systems 
that demonstrate spin polarization from unpolarized in- 
put, in many cases the effect is subtle J^ii^ii^!^ Here we 
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give an idealized thought experiment showing that an un- 
polarized input current can produce a spin polarized out- 
put current. Consider a system with = 1 and M = 2, 
as illustrated in Fig.[T] All input electrons are incident on 
a perfect skew scatterer which sends spins quantized in 
the +z direction into exit channel 1 and spins quantized 
in the — z direction into exit channel 2. Exit channel 2 
has a region with Rashba spin-orbit coupling^ which is 
precisely of the strength and length necessary to rotate 
—z spins to +z. Thus, all spins incident from the left lead 
exit with their spins up, and the system respects TRS, 
since skew scattering and Rashba spin-orbit interaction 
are each time reversal symmetric. 

We illustrate by constructing S explicitly. We can ex- 
press the scattering matrix for this thought experiment 
(up to an overall phase) in the 6x6 and 3x3 represen- 
tations as 



Using Eq. m the spin conductance is 
In particular, if has no quaternion part, then = 
0. We have 



out _ 1 - SkkS 
Wkk- ^ 



'KK 



(15) 



and SkkS*j^j^ is real, so € M and = 0. This proof 
applies with channels that are fully open or have tunnel 
barriers, as it requires only that the S-matrix satisfy TRS 
and unitarity, which are unchanged by tunnel barriers. 

We note further that ii K > 2 then 1) the reflected 
current in any of the K — \ input channels can be spin 
polarized, and 2) if the input current goes through less 
than K — 1 channels, then the remaining channels can 
have a spin polarization, as shown in the example of Fig. 

m 



5 = 








(10) 

(11) 




where 9 e [0, 27r) and r and t have been determined by 
the above description, while the rest of the matrix is given 
by TRS and unitarity. The unpolarized input quaternion 

density matrix is = ( ^ o \ giving 



(12) 



so Eq. 2] gives ~ p — z, as stated above. 

We now prove that having at least two channels in 
the outgoing lead is essential. That is, for a dot with 
TRS and K channels in attached leads, if an unpo- 
larized equally weighted incoherent current is sent into 
N = K — 1 oi the channels, then the spin polarization 
in the remaining channel must be zero. This result has 
been shown beforCf^ but the quaternion formalism with 
density matrices makes it particularly transparent, so we 
include the proof here. 

We start with 



w 



J_ (In 
2N \ 



1 



K 



K 



2N 



(13) 



where Pk is the projection onto the K^^ channel. The 
quaternion scattering matrix satisfies S — since TRS 
is unbroken, and 



SS'f - SPkS'' Ik - SPkS'' 



2N 



2N 



(14) 



Note that S = implies both = S* and S^^ £ C for 
i = l...K. 



III. RANDOM MATRIX THEORY 

We estimate the expected spin polarization in realistic 
situations by using random scattering matrix theory. We 
assume that the mean dwell time of particles in the dot 
is much greater than the Heisenberg time th = 27r?i/A, 
where A = 2Trh^ /mA is the mean orbital level spacing, 
m is the effective mass, and A is the area of the dot. 
We further assume the strong spin-orbit limit, where the 
spin-orbit time Tgo is much less than r^. For a chaotic 
quantum dot, Td = mA/HK, where K is the number 
of fully open orbital channels attached to the dot, so 
for a sufficiently large A, even a material with "weak" 
spin-orbit coupling will be in the strong spin-orbit limit. 
The crossover from weak to strong spin-orbit coupling in 
chaotic quantum dots has been studied in the K ^ 1 
limit in the context of adiabatic spin pumping.^^ 

For dots with strong spin-orbit coupling, we assume 
that the S-matrix is chosen from the uniform distribu- 
tion of unitary matrices subject to TRS, called the cir- 
cular symplectic ensemble (CSE)i ^^i^° We find the root 
mean square (rms) magnitude of the spin conductance 
on averaging over the CSE, which gives the typical spin 
conductance magnitude to be expected from chaotic de- 
vices. Such an averaging can be realized in practice by 
small alterations of the dot shape.— 

By symmetry, (g*) = for fi — 1,2,3. Using Eq. IH 
we evaluate 

((/)') = 4(tr(a,,Pfl5u;'"5t)tr(a^Pfl5«;'"5t)) ^ (ig) 

where we sum over /i. 

We use the technique for averaging over the CSE 
described by Brouwer and Beenakker in Section V of 
Ref. [s^. We need just two generic averages, which we 
will use repeatedly. The first is of the form {Fi{S)) = 
{tT{ASBS^)') , where A and B are constant KxK quater- 
nion matrices and the average is taken over S chosen from 
the CSE oi K X K quaternion self-dual matrices. Then^^ 

^^'^ = 2^[2tr(A)tr(i?) - tr(Ai?«)]. (17) 



4 



The second average we need is {F2{S)) = 

{tT{ASBS'')tiiCSDS'')) where A, B, C, D 

are constant K x K quaternion matrices and 
AB = AD = CB = CD = 0. We find^^ 

(F2) = - l)[4trAtrBtrCtr7:> + tr(ylC)tr(Bi:i)] 

- [trAtrCtr(BD) + tr(AC)trmri:>] } , (18) 



where A = K{2K - 1){2K - 3). 
Using Eq. [THl we find 



^M{M - 1) 



(19) 



where we used tr(CT^Pfl) = 0, tr(P|) = M, trw'" = 1/2, 
and tr[(u;'")2] = 1/4A^. Note that when M = 1, 
{{d^y} = consistent with the general symmetry. 

If we are interested in the mean square polarization of 
the exit current, (p^) = {{g^)"^ g~'^) , we can approximate 
it by ((g'')^) / (g)^- This approximate form is useful for 
analytical progress and will be compared to numerical 
results. Using Eq. [T71 



(5) 



2M 



2K 

which, combined with Equation 1191 gives 



3{M - l)i2K - 1) 
UINK{2K-3) ' 



(20) 



(21) 



We study 

\2 



the approximation ((5*) ^5 ^) ~ 

((5'*)^) / (5)'^ numerically. We choose a 2K x 2K 
complex Hermitian matrix from the Gaussian unitary 
ensemble^^ and find the unitary matrix U which di- 
agonalizes it. We multiply columns of U by random 
phases, map U into a K x K matrix of quaternions, 
and construct unitary self-dual S by setting S = UU^, 
giving 5* chosen from the CSEi^ 

Figure [5] shows the numerical and analytical results, 
which agree quantitatively for {g^} and (((7^)^) and qual- 
itatively for (p^)- The largest percentage disagreement 
for (p2) is 7%. 




FIG. 2: (color online) Numerical (symbols) and analytical 
(lines) results for normalized mean conductance {g), rms spin 
conductance g", and rms spin polarization p of current exiting 
a chaotic quantum dot with A'^ (M) channels in the entrance 
(exit) lead. An average over 60000 S-matrices from the CSE 
was performed for each data point. The lines are from Eqs. 



the charge/spin current evenly between the channels, 
and thus models dephasing processes that preserve elec- 
tron spin. First, consider 77^ = tT{af^Pci,Sw™S^), where 
/i = 0, 1, 2, 3, CTo = I2, Pci> is the projection operator onto 
the dephasing lead's channels, and is the input den- 
sity matrix. Then 2r]Q is the probability for a particle to 
enter the dephasing lead, and 2ff is the spin conductance 
into the dephasing lead, which is proportional to the spin 
current into the dephasing lead. 

We reinject from the dephasing lead with 



IV. DEPHASING 

We add dephasing to this setup using the dephasing 
voltage probe technique.— 1^1^ We add a fictitious volt- 
age probe drawing no current with = 2T:h/ At^ fully 
open orbital channels, where is the dephasing time. 
We extend this model to preserve the spin of the rein- 
jected electrons. 

In contrast with previous work, we explicitly model 
reinjection of electrons from the voltage lead by modify- 
ing to include incoherent reinjection from the dephas- 
ing lead. The reinjection matches the total charge/spin 
current absorbed by the dephasing lead, but distributes 



'■^C' -i"'''" <^^) 



where we sum over repeated index ^. We set = VtJ./^'P^ 
which ensures the reinjected charge/spin current equals 
the absorbed charge/spin current. Some of this reinjected 
current reflects back into the dephasing lead, so it must 
be reinjected again. We define a 4 x 4 complex matrix 
0/^1/ — ir{a^Pci,Sa^P^S''), which gives the charge/spin 
current in the dephasing lead due to this reinjection. 
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FIG. 3: (color online) Numerical (symbols) and analytical 
(lines) results for normalized mean conductance {g), rms spin 
conductance g", and rms spin polarization p of current exiting 
a chaotic quantum dot with A'^ = 1 channel in the entrance 
lead. M and A''^ are the numbers of channels in the exit 
and dephasing leads, respectively. An average over 60000 S- 
matrices from the CSE was performed for each data point. 
The hues are from Eqs. [27] and [28l 



Defining w'" = w^^ + w™, this procedure gives 

C30 



(6 



n-l\ 



= P0(T^tr(a,P^5<5^)(iV^(5^, - e^,)-\ (23) 

wiiere we sum over repeated indices ^^v — 1, 2, 3. This 
result holds for any input current, not just the unpolar- 
ized incoherent Wq" discussed here. 

We approximate w^^ by replacing with its average 
in Eq. [11 similar to Eq. [HI Using Eq. [T71 

(0M-> = ^ , [2(iV0 " l)<5^o^.o + 5,..], (24) 

where = N + M + N^. We further replace 

ir[afj_Prj,Sw™S'^) by its average, 

(tr(a^P^5«;-5t)) = 4^^, (25) 



2i^^ - 1 



which gives 



This turns out to be the same result as if we had chosen 
= in Eq. [Hjfor = 1, 2, 3. That is, in the approxi- 
mation of Eqs. [24H251 if we have total spin decay in the 
dephasing lead, then Eq.[2niis unchanged. Note that Eq. 
[26] satisfies unitarity only on average; the total probabil- 
ity of exiting either through the right or left lead equals 
1 only on average. 

In this approximation, we use Eq. [TS]to find 



where = K42K^ - 1){2K^ ~ 3). 

Note that even if there is only one outgoing channel, 
M — 1, the spin conductance is predicted to be nonzero 
due to dephasing. The dephasing induced spin conduc- 
tance is present for > 1, as shown in numerical sim- 
ulations in Fig. [3] and is not an artifact of Eq. [^S] In 
the case of = 7\/ = 1 , an exact treatment shows that 
— 0, contrary to Eq. 1271 even with arbitrary tunnel 
barriers between the leads and the sample. As shown in 
Fig. [21 Eq. [27] works well for M > 1 or A^^ > 1. 

We can modify this model to have Ncj, dephasing leads 
each with one channel, each separately reinjecting the 
same charge/spin that it absorbs. In this model, too, a 
nonzero g'^ can be produced for > 1 (not shown). 

Brouwer and Beenakker modified the third-lead de- 
phasing model to make dephasing uniform in phase space 
by placing a tunnel barrier with transparency F between 
the dephasing lead and the dot, with F ^ and A^^ oo 
while maintaining FTV^ = 2Trh/ At^.^^ The S-matrix is 
then not drawn from the CSE, and simple analytical 
results in the spin-orbit coupled system are challeng- 
ing. We study this model numerically and find that for 
fixed r^, it gives qualitatively similar results to the sim- 
pler model described above; in particular it also gives a 
nonzero spin current when M = 1 (not shown). Without 
a microscopic model of dephasing, it is possible that this 
dephasing induced spin current with Af = 1 is an artifact 
of third lead dephasing models, but all three variants of 
third lead dephasing discussed here show this effect, so 
dephasing gives a loophole for producing spin currents 
even when M = 1. 

Returning to the single dephasing lead with F = 1 , we 
estimate {p^} as above, where we modify (g) to include 
the dephasing lead. Using Eqs.fTTlandl^ this gives 



(.9) 



2MK,, 



K{2K^ - 1) ' 



(28) 



w'" « 4" + P^/2K. 



(26) 



We estimate {p^) « {{g'f) / {gf, using Eqs. [17] and [M] 
Comparison of these approximations to numerical eval- 
uations is shown in Fig. [3] Again we find that the nu- 
merical and analytical results agree qualitatively, except 
when N^ = M = 1. 
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V. FINITE TEMPERATURE 

If the temperature T > A, the polarization wiU be 
further suppressed by electrons of different energy feel- 
ing uncorrelated scattering matrices. This effectively in- 
creases the number of orbital channels, which decreases 
the residual polarization. We consider unpolarized inco- 
herent flux from the left lead at temperature T . Adapting 

Datta,.2a we take w'"(e) = -§7 w ( ^" Om ) ' ^^^^^ ^^'^^ 
is the Fermi distribution. If the scattering matrix for par- 
ticles of energy e is S'(e), then w™*(e) = S'(e)w'"(e)5't(e). 
We approximate S{e) as correlated only within energy 
intervals of scale A (see Ref. HO for an equivalent treat- 
ment). That is, we take 

{Sab{e)Sl{e')) = AJ(e - e') {Sab{e)sUe)) , (29) 

and 

{Sab{e)SUe')Slj{e)Sl^{e')) 

= {Sa,{e)Slf{e)) {s.,{e')Sl^{e')) (30) 
+ A<5(e - e') {Sa,{e)SUe)Slf{e)Sl^{e)) , 



which are valid only for T ^ A, which is often true for 
chaotic quantum dots. For T « A, (5(6)5^(6)) can be 
calculated using the random Hamiltonian method4i 



We need an average over a new function. 



h{e,e') = ne)ne')iv[AS{e)BS\e)]iv[CS{e')DS\e% 



where AB ^ AD = CB = CD = and /' = df/de. We 
evaluate the average of h with the K x K quaternion ma- 
trix S{e) chosen from the CSE along with Eq.[3Dl giving. 



I dede' {h{e,e')) = ^--l-^tr(A)tr(i?)tr(C7)tr(I?) 
A f 

+ — / def'{ef{{K - l)[4trAtrmrCtri:> + tr(AC)tr(SL>)] - tryltrCtr(B£i) - tr(ylC)trBtri:>} 
" (2jrll)2 tK^)tr(i?)tr(C)tr(J) 
+ ^-{{K - l)[4tryltrmrCtrD + tr(^C)tr(BD)] - trAtrCtr(Bi:)) - tr(AC)trBtri:>}. (31) 



Using Eq. [3T]in place of Eq.[T8l we evaluate ((5*)^) as 
above, which simply multiphes Eq. [T9]by Also, (g) 
is unaffected by temperature, so Eq.[2I]is also multipUed 
by A/6T. 

When dephasing and temperature are both included, 
the scattering matrix is correlated on the scale of the 
level broadening, A(l -I- so we replace A in Eq. 

[Sniby A(l + N^/2). Eq. [^is then multiplied by A(l + 
N^/2)/6T, and Eq.[15]is unchanged. 



VI. DISCUSSION 

This spin polarization should be able to be produced 
and detected experimentally. Even quantum dots in n- 
type GaAs/AlGaAs heterostructures have been observed 
to have sufficiently strong spin-orbit coupling to approach 



the RMT symplectic limiti^'^ If the spin-orbit coupling 
is not strong enough for the S-matrices of the dot to 
be drawn from the CSE, the spin polarization predicted 
here will be reduced but should still be present. In a 
given material with fixed spin-orbit coupling strength, a 
sufficiently large quantum dot will be well described by 
the CSE, with a possible increase in dephasing rate as 
the dot size increases. 

At zero temperature, our discussion has assumed that 
all electrons passing through the dot see the same S- 
matrix, which is valid when the applied potential dif- 
ference is less than the mean level spacing A. Since 
A cx A"^, as the dot area is increased to approach the 
strong spin-orbit limit, the window of voltages where 
these results apply shrinks. The effects predicted in 
this paper are most likely to be observable in a material 
with inherently strong spin-orbit coupling, such as p-type 
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III/V hctcrostructurcs. Note that ((,g*)^) = ((.9")^) /3 
and (p^) = (p^) /3, so if a measurement technique or 
application is only sensitive to spin polarization along 
a particular axis, then the rms predictions for the ^- 
component of the polarization and spin conductance are 
only VS times smaller than the results stated above. 

We have shown that quantum dots with spin-orbit cou- 
pling can generate spin polarized currents without mag- 
netic fields or ferromagnets, except with only one out- 
going channel and TRS, when such a device cannot pro- 
duce a spin current. Mesoscopic fluctuations can be large 
enough to give appreciable spin currents in devices with 
a small number of propagating channels. Even if the 
spin-orbit coupling is weak, a suSiciently large device will 



show these eff'ects. 
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